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ABSTRACT 

A phase transition occurring in the inner core of a neutron star could be as- 
sociated to a density discontinuity that would affect the frequency spectrum 
of the non-radial oscillation modes in two ways. Firstly, it would produce a 
softening of the equation of state, leading to more compact equilibrium config- 
urations and changing the frequency of the fundamental and pressure modes of 
the neutron star. Secondly, a new non-zero frequency g- mode would appear, 
associated to each discontinuity. These discontinuity g-modes have typical 
frequencies larger than those of g-modes previously studied in the literature 
(thermal, core g- modes, or g-modes due to chemical inhomogeneities in the 
outer layers), and smaller than that of the fundamental mode; therefore they 
should be distinguishable from the other modes of non radial oscillation. In 
this paper we investigate how high density discontinuities change the fre- 
quency spectrum of the non-radial oscillations, in the framework of the gen- 
eral relativistic theory of stellar perturbations. Our purpose is to understand 
whether a gravitational signal, emitted at the frequencies of the quasi normal 
modes, may give some clear information on the equation of state of the neu- 
tron star and, in particular, on the parameters that characterize the density 
discontinuity. We discuss some astrophysical processes that may be associated 
to the excitation of these modes, and estimate how much gravitational energy 
should the modes convey to produce a signal detectable by high frequency 
gravitational detectors. 
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2 G. Miniutti et al. 
1 INTRODUCTION 

The equation of state (EOS) of dense matter at supranuclear densities has been a long 
sought-after problem in nuclear physics. Since the first observational evidences of their ex- 
istence in the late 60's, neutron stars (NS) have proved to be a unique astrophysical envi- 
ronment in which very different fields of physics, ranging from nuclear physics to particle 
physics and general relativity, can be tested. Accurate estimates of some properties of NSs, 
such as masses and radii, associated to further information on the thermal history, would 
allow to set constraints on the nature of nuclear forces. This has been attempted by several 
means (but so far with limited success), by correlating mult i- wavelength observations with 
different telescopes, X-ray, or 7— ray observatories ( Thorsett fc Chakrabarty 1999| ; [Rutledge 



\et al. 200T| ; Pons et al. 2002E| ), or through the study of neutrinos from proto- neutron stars 



( Pons et al. 1999) ; Pons et al. 2001| ), or studying QPO properties QOsherovich fc Titarchuk 



[I999t [Stella fc Vietri 1999Q . 



In the next future, a new observational window will be opened by a network of gravita- 
tional wave detectors: the ground based interferometers (LIGO, VIRGO, GEO600, TAMA), 
which are now in the final stage of construction or in the commissioning phase, will explore 
the frequency region of 10-1000 Hz, and the space-based interferometer LISA will enlarge 
the observational window to 10" 4 — 10 _1 Hz. In addition, a joint panel of European experts 
have performed an assessment study for a new generation of interferometric detectors with 
extremely high sensitivity in the kHz region (EURO, EURO-XYLOPHONE), which would 
provide a very powerful mean to investigate the physics of neutron stars through gravita- 
tional waves. Since these instruments are currently under study, it is important to envisage 
which sources they would be able to see, and what kind of information could be inferred from 
the detected signals. In this paper, we address to one interesting possibility in this exciting 
new field: the prospects of probing phase transitions occurring in the inner core of neutron 
stars. In particular, we shall focus on phase transitions that produce a density discontinuity. 

Density discontinuities may occur either at high or low density, depending on the physi- 
cal processes prevailing in the different regions of the star. At low density, the characteristic 
chemical profile in the outer layers is determined by the history of the star, because shell 
burning, flash nuclear burning and accretion phenomena leave layers of different composition 
on its surface. At the interfaces between different layers, the chemical composition changes 
abruptly, and if no significant diffusion is present, the density gradient is well approximated 
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Discontinuity modes in NS 3 

by a discontinuity. This introduces an additional local source of buoyancy which gives rise to 
a discontinuity g-mode, even if no other source of buoyancy is present, i.e. even in zero tem- 
perature NSs. The properties of the g-modes due to discontinuities in the outer layers in cold 
NSs have been studied by Finn (1987) and McDermott (1990). Subsequently, Strohmayer 
(1993) computed the oscillation frequencies of two finite-temperature models with a density 
discontinuity associated to the 56 Fe to 62 Ni transition. He found that if the discontinuity is 
trapped into the crust no mode can be found directly associated with it, while if it is located 
in the fluid ocean a g-mode appears, and the frequency distribution of the g-modes due to 
entropy gradients changes. 

However, density discontinuities are not necessarily confined to the outer layers of a neu- 
tron star: they could also arise as a consequence of phase transitions in the inner core, where 
the equation of state is still poorly known. Several first or second order phase transitions have 
been proposed to occur at supranuclear densities, and they typically involve pion and/or 
kaon condensation as well as the transition from ordinary nuclear matter to quark matter 
( [Hciselberg fc Hjorth- Jensen 2000| [Prakash et at. 1997 and references thereTri| ). Whether 
these transitions are first or second order, and whether they admit density discontinuities 
or a continuous transition with a mixed phase, is still matter of debate. In modeling the 
phase transition, one possibility is to impose local conservation laws and to use a Maxwell 
construction, which produces a density discontinuity. Glendenning (1992) pointed out that, 
since conservation laws must be imposed globally, any first order phase transition should be 
accompanied by a smooth mixed phase region, which satisfies the Gibbs rules. However, it 
should be noted that there is much uncertainty on the effect that Coulomb repulsive forces 
and surface tension have on the structure of the mixed phase. If Coulomb and surface ener- 
gies are large, the mixed phase is disfavoured, the region where it occurs shrinks, and this 
seems to be the case in a quark-hadron phase transition ([Afford et at. 200T ). In this case, 
density discontinuities may appear near the boundaries of the Gibbs' structured phase, when 
the volume fraction of one of the phases is very small. 

Astrophysical observations play an important role in clarifying these issues, because 
neutron stars provide a unique opportunity to study the behaviour of matter in this very high 
density regime. In the present paper we study how density discontinuities affect the spectrum 
of the quasi normal modes of a neutron star, with the purpose of understanding what kind 
of information can be inferred on the high density equation of state from the gravitational 
signal emitted by an oscillating neutron star. We shall consider some astrophysical processes 
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4 G. Miniutti et al. 

in which the quasi normal modes could be excited (a glitch, the onset of the phase transition, 
binary coalescence of NS's) and we shall estimate the amount of energy that should go into 
such modes to produce a signal that may be detected by the proposed high frequency 
gravitational detectors. In this respect, our work differs from a previous work of Sotani et 
al. ( [Sotani, Tominaga fc Maeda 200"T| ), who focused on the calculation of the modes and 
discussed the stability of the considered models. For simplicity, we consider the case of cold 
neutron stars and assume a simple polytropic EOS. 

In Sec. 2 we briefly describe the relevant equations of stellar perturbations which we 
integrate to find the mode frequencies. In Sec. 3 we discuss the results of the integration, 
extracting the information that the discontinuity g-modes carry on the parameters of the 
density discontinuity. In Sec. 4 we consider some astrophysical processes that are associated 
to the excitation of the quasi-normal modes, and discuss the detectability of the emitted 
gravitational radiation. In Sec. 5 we draw the conclusions. 



2 THE MATHEMATICAL FRAMEWORK 

In order to find the frequencies of the quasi-normal modes, we will integrate the equations 
describing the polar, non radial perturbations of a non rotating star as formulated by Lind- 
blom and Detweiler ( Lindblom fc Dctwcilcr 1983 ; Dctwcilcr fc Lindblom 1985 ). We write 



the perturbed metric tensor as 

ds 2 = e v {\ + r l H 0lm Y lm e^)dt 2 - 2iuor l+l H llm Y lm e luJt dtdr + e A (l - r l H 0lm Yi m e iuJt )dr 2 

+ r 2 (l-r l K lm Y lm e i " t )(d6 2 + sin 2 ed<j) 2 ), (1) 

and the polar components of the Lagrangian displacement of the perturbed fluid elements 

as 

f = e- x / 2 r l - l W lm Y lm e^ , 

i e = -r l - 2 V lm d e Y lm e^ , (2) 
= -r l (r sin ey'V^Y^ . 

Yi m (6, <j)) are the spherical harmonic functions. In the following we shall restrict our analysis 
to the I = m = 2 component, which dominates the emission of gravitational radiation. By 
defining the variable X 

X = u 2 (p + p)e~ u/2 V - r- l e^- x)/2 p'W + \{p + p)e u/2 H , (3) 
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Discontinuity modes in NS 5 

where p and p are the fluid's pressure and energy density and the prime indicates differ- 
entiation with respect to r, one can write a fourth-order system of linear equations for 
(H h K,W,X) 

H[ = ^(-[(l + l) + e\2M(r)/r + r 2 (p-p))]H 1 + e x [Ho + K-4(p + p)V] 
K' = 1 (H + {n+ 1)^i- [(I + 1) -r^u']K -2(p + p)e x/2 W 



W 

X' 



(1+1) 



x/2fe~ u/2 X l{l + l)V 1 \ 
W + re A/2 | i ^— + -H + K ) , 



IP 



l -X + (p + p)e^ {I (i - i „') # + l[rV + (n + 1)]^ + ri/' - 1)# 
,Z(Z + 1)V 1 



e A/2 [(p + p) + Cj 2 } - l -r\r- V^V) ' 



W 



(4) 



where n = (Z — 1)(Z + 2)/2, = w e and 7 is the adiabatic index 
p + p (dp\ P + P 2 



p \dp 



P 



(5) 



The functions V and Ho are 



2r 2 e -^/2 X _ [I ( n + l) rzy ' _ r 2 cD 2 ] e - A iJi + [n - w 2 r 2 - § u'(3M(r) - r + r 3 p)]K 



3M(r)/r + n + r 2 p 



J 2 V(P + P) 2r 2 °. 



(6) 



For each given value of Z and u, Eqs. (^) admit only two linearly independent solutions 
regular at r = 0; we find the general solution as a linear combination of the two, such that 
the perturbation of the Lagrangian pressure, Ap, vanishes at the surface. This procedure is 
slightly different from that proposed by Lindblom and Detweiler (1983), but a comparison 
of the results of the two procedures shows that they are absolutely equivalent. 

It should be noted that, in this formulation, both the £ r and Ap are continuous even 
when the density is discontinuous. 

It is well known that outside the star, the variables associated to the fluid motion are 
zero and the equations reduce to a second-order equation for the Zerilli function Z 
( d 2 



dr'i 



UZ , 



u 



2(r - 2M) 



n 2 (n + l)r 6 + 3MnV + 9M 2 (M + nr) 



(7) 



r 4 (W + 3M) 2 

where r* = r + 2Mlog(r/2M — 1) and M is the mass of the star. The asymptotic behaviour 
of the Zerilli function is 



A in (u)e iuJI * + A out (u)e- 



(8) 
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A quasi normal mode of the star is defined to be a solution of the perturbed equations 
belonging to a complex eigenfrequency uj q + icj im , which is regular at the center, continuous 
at the surface, and which behaves as a pure outgoing wave at infinity. Therefore, if ujq is 
the mode frequency, A in (uo) must vanish. In order to find the frequency of the quasi-normal 
modes we use the following algorithm flUhandrashekar fc Ferrari 1990| ), which is appropriate 
when the real part of the frequency ujq is much larger than the imaginary part u)\ m = 1/r, as 
it is in our case. For real values of the frequency, the function Z must have the asymptotic 
form 



a 



n + lp 



uj r 
n + 1 a 
uj r 



2uj 2 
1 

2^ 



1 

^7 + 



cos o;r* 



a 



+ 



sin cur* 



(9) 



n(n + l)a - 3Muj [1 + ~)P 

n(n + 1)0 + 3Muj (l + 

The functions a(u) and (3{uj) can be determined by matching the numerically integrated 
solution with the above asymptotic expression for Z. The amplitude of the standing grav- 
itational wave at infinity is (a 2 + /3 2 ) 1//2 , and it can be shown to have a deep minimum at 
the frequency ujq where A^iujo) vanishes ( |Chandrasekhar, Ferrari fc Winston 1991| ); thus, we 
find the frequency of a quasi normal mode by searching for the values of ujq where (a 2 -!-/? 2 ) 1 / 2 
has a minimum. 



3 RESULTS 

To study the effect of high density discontinuities on the oscillation spectrum of a NS, we 
shall consider a simple polytropic equation of state of the form 

Kp r p > p d + Ap 

P=\ ( ApV r > (10) 

K{l + -^-J p r p<p d 

where the discontinuity of amplitude Ap, is located at a density pd, which we choose to be 
greater than the nuclear saturation density po = 2.8 ■ 10 14 g/cm 3 . This simplified equation of 
state will allow us to study in a systematic way the properties of the quasi normal modes, 
and the relation between the mode frequencies and the global properties of the chosen 
models. This EOS has been used in the past to study the g-modes associated to low density 
discontinuities ( |Finn 19871 [McDermott 1990] ) and, in a slightly modified form, to study the 



g-modes due to high-density discontinuities (Sotani et al, 2001). 

To find the frequency of the surface g-modes, Finn introduced the so called slow-motion 
formalism ([Finn T986| ) which is very accurate in the low- frequency limit. In our case the 
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Discontinuity modes in NS 7 

mode frequencies are higher since they are associated to high density discontinuities, and 
the full system of perturbation equations has to be integrated. As a consistency check of 
our code, we have repeated some of Finn's calculations for low density discontinuities ( |Finn 
1987] ) finding an excellent agreement (better than three significant digits). 



It should be mentioned that g-modes may have real, imaginary or zero frequency, de- 
pending on whether the considered stellar model is stable under convection; different regimes 
can be identified by looking at the sign of the square of the Brunt- Vaisala frequency 

JV^Ie-VLfl-i) , with 7o = ^±pil. 

2 V \1 7o/ VP 

A real, imaginary or zero N correspond, respectively, to convective stability, instability or 

marginal instability. If the star is isentropic and chemically homogeneous, 7 = 70 and all 

g-modes degenerate to zero frequency. Thus, in the case of our zero temperature, chemically 

homogeneous star, N is zero everywhere, except when the density changes abruptly at 

r = Rd, at that point, the adiabatic index 7 remains finite, the equilibrium index 70 vanishes, 

N is different from zero, and a g-mode appears. 

Before discussing the behaviour of the quasi normal mode frequencies, let us shortly 

describe how the structure of a neutron star changes if a phase transition associated to a 

density discontinuity occurs in its inner core. In Figure [I] we plot the mass-radius relation 

for a star with and without discontinuity. The parameters are Pd = 7 • 10 14 g/cm 3 and 

Ap/p d = 0.3. The polytropic index is T = 2 in both cases, and we set K(l + Ap/p d ) 2 = 

180 km 2 , in order to have the same equation of state for the two models when p < p d . Figure 

D shows that the maximum mass is lower for the model with a density discontinuity; this is 

a general behaviour, which is due to the softening of the equation of state introduced by the 

discontinuity. As a rule of thumb, for a fixed mass, a star with a discontinuity is more compact 

than one with the same mass and no discontinuity at all; for instance, from Figure f| we see 

that for M = 1.4 M the star with no phase transition has a radius of R = 13.44 km, while 

the one with a density discontinuity has R = 9.66 km. These general properties continue 

to hold also if EOSs more realistic than eq. (|T0|) are used, and if phase transitions due to 

kaon/pion condensation or to the nucleation of quark matter are considered ( [Heiselberg & 

[Hjorth-Jcnsen 20"0T ). Thus the EOS we use, though approximate, allows us to infer the mode 

properties that are known to depend much more on the star's macroscopic properties (mass, 

radius, etc.), than on the specific microscopic interactions. 
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8 G. Miniutti et al. 

In our study, we shall consider stars with a fixed mass M = 1.4 M , because astronomical 
observations show that NSs that have not spent a significant part of their life accreting matter 
from a companion in a binary system, have masses in a narrow range around this fiducial 
value ( frhorsett fc Chakrabarty 1999| ). We shall choose different values of the polytropic 
exponent T = (1.67, 1.83, 2.00, 2.25), and assign the constant K in such a way that models 
with the same V have the same EOS for p < pd , independently of the particular values of 
Ap/pd and pd- Since we require that all stars have the same mass, we will have to adjust the 
central density accordingly. This means that different stars have a different matter content 
in their core. Furthermore, we select models with reasonable radii (8.5 km < R < 17.9 km) 
and, since we want to study high density discontinuities, we fix their location at < 0.8 R. 
The amplitude of the discontinuity is chosen to vary within Ap/pd = (0.1 — 0.3); larger 
values would not allow stable stellar models with reasonable values of mass and radius. 

In Tables |l]-[| we summarize the results of the numerical integration of the Eqs. @ and 
(0) for these stellar models, as well as the equilibrium parameters. The frequencies of the 
fundamental mode (i/f), the first pressure mode (z/ p ), and the discontinuity g-mode (f g ) are 
tabulated for different values of the parameters. The models considered in Table [I] have 
the same value of the polytropic exponent, T = 2, and refer to three different values of 
Ap/pd = ( 0.1,0.2,0.3 ). In the remaining three Tables (§-§), we change the polytropic 
index and for each value we consider a different Ap/pd. These choices allow to study the 
mode frequencies for a wide range of stellar models. 

The dependency of the mode frequencies on the variables relevant to our analysis is better 
clarified in the following figures. In Figure |2], we plot the frequency of the fundamental mode 
and of the first p-mode as a function of pd, for models with T = 2 (Table |I]), and for the 
three selected values of Ap/p d . We choose the T = 2 models because they are sufficiently 
representative of the general behaviour. The curves representing z/f and u p for different 
Ap/pd tend to a limiting value (shown as a black dot in the Figure), which is the frequency 
of the modes when there is no discontinuity inside the star, i.e. when pd exceeds the central 
density. Both u { and z/ p increase with Ap/p d , and decrease with p d . This behaviour is 
not surprising. Indeed, from the data in Tables we see that, having fixed the mass of 
the NS, both the mean density and the compactness increase with Ap/pd and decrease 
with pd; from the Newtonian theory of stellar perturbations it is known that the f-mode 
frequency is proportional to JM/R 3 , and the first p-mode frequency scales with the star's 
compactness Mj R; thus, the results shown in Figure 2 confirm this general behaviour, despite 
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Discontinuity modes in NS 9 

the structural changes produced by the density discontinuity. The dependency of u { and 
Up on the mean density and compactness of the star has been studied in general relativity 
for a large number of realistic EOSs; empirical relations have been derived between the 
mode frequencies and the parameters of the star, that could be used to infer the mass and 
the radius of neutron stars if a gravitational signal emitted by these pulsation modes were 
detected ( |Andersson fc Kokkotas 1996|; |Andersson fc Kokkotas 1998| ; [Kokkotas, Apostolatos 
Andersson 200l|) . These fits were obtained on the assumption that the density is continuous 
inside the star, and it is interesting to see whether a discontinuity introduces deviations from 
the expected behaviour. For instance, in Figure 3 we plot, as a continuous line, the linear fit 



which relates and yM/i? 3 obtained by Andersson and Kokkotas (1998) (AK98-fit). In 
the same figure, we also plot the values of z/f which we find for models with different values 
of r and lS.pl pa- The frequencies of the models with Ap/ p& = are shown as black dots. We 
see that if a density discontinuity occurs in the inner core, V[ deviates from the fit; the fit 
is expected to have an error less than 1% ( |Kokkotas, Apostolatos fc Andersson 2001| ), but 
the deviation introduced by a density discontinuity is certainly larger. Thus, the detection 
of an f-mode frequency would still allow to estimate the mean density of the star, but with 
a larger uncertainty. Similar conclusions can be drawn by plotting the p-mode frequency as 
a function of Mj R. 

Let us now focus on the properties of the g-mode which is produced by the density 
discontinuity. The possible values of u g cover a wide range, v g ~ [ 0.5 — 1.4 ] kHz (see Tables 
and they are the highest frequency g-modes that one can expect in neutron stars. 
Indeed, g-modes due to finite temperature effects have very low frequencies ( pVIcDermott, van 
Horn & Scholl 19831 IMcDermott, van Horn fc Hansen 1988]) , while those due to composition 
gradients in neutron star cores ( [Reisenegger fc Goldreich 1992| ; [Lai 1994]) and to density 
discontinuities in the outer layers ( Finn 1987 ; McDcrmott 1990 ) have, typically, v < 200 Hz. 
Thus, should a g-mode be detected, one could establish if it is due to a discontinuity at 
supranuclear densities by looking at its frequency. The observation of such a g-mode would 
also carry information on the parameters that characterize the discontinuity. This is clear 
from Figure |], where we plot z/ g as a function of pd, for all the considered models. We see 
that the g-mode frequencies belonging to the same density jump, Ap/ pd, cluster in a definite 
region of frequency, independently of the other parameters of the EOS (r and pd), and this 
is a remarkable feature. 
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10 G. Miniutti et al. 

Furthermore, for a fixed value of Ap/p^, v z exhibits a maximum as a function of the 
location of the discontinuity. To better show this behaviour, in Figure |5] we explicitly plot 
u g , computed for two stars with the same Ap/p^ = 0.2 and having T = 1.83 and T = 2, 
respectively, as a function of the radius at which the discontinuity occurs, R&. The reason why 
there is such maximum is the following. In a chemically homogeneous, zero temperature star 
with no density discontinuity, the g-mode spectrum is degenerate at zero frequency. When 
we consider stars with very close to center, or when we put the discontinuity very close 
to the surface, we are approaching the model with no discontinuity, and in both cases u g 
tends to zero. Thus a maximum has to be expected at some Rd , intermediate between these 
two extreme cases. For all the models we studied the maximum g-mode frequency is found 
when 0.4 < R d /R < 0.6. 

As already shown by Finn (1987) following a Newtonian analysis (|Landau fc LifshTti 
|1959|) , if the discontinuity is in the low-density region near the surface of the star, the 
dispersion relation for the oscillation modes can be written as 

where p = M/R? is the mean density of the star, Ar = R — R^, and I is the order of 
the considered multipole. This relation is obtained by solving a stratified two-fluid problem, 
taking the lower density fluid to be infinitely deep and approximating the spherical symmetry 
of the star by plane symmetry, assumptions that are reasonable if the density discontinuity 
is near the surface of the star. Obviously, they are no longer satisfied when the discontinuity 
occurs at high density and R d ~ 0.5 R. However, we find that the maximum g-mode 
frequency for each assigned EOS satisfies a similar relation 



u™ x 2 ~/(/ + l)C 2 ^ , (13) 
where 



which differs from the corresponding term in eq. (|T2"|) only because the mean density is 
replaced by the mean density of the "internal" core (p > pa), Pi- Considering all our stellar 
models, we find the following linear fit 

v™ x = \ 0.070 + 0.905 • ( T ) 1 kHz (15) 

g L \0.025 knrV J V ; 

with an error smaller than 2%. For our models, 0.016 km -1 < C < 0.037 km -1 . This linear 
behaviour has an interesting consequence. Suppose that a g-mode with a given frequency 
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Discontinuity modes in NS 11 

z/ obs could be detected. Then, from Figure [| we would be able to constrain the amplitude of 
the density discontinuity Ap/p^, and from the fit ( |I5|) we could set a lower limit on C and 
consequently on p\. All the EOSs with C < C{v hs) would be ruled out because they would 
have a g-mode with a maximum frequency lower than the observed one. 



4 EXCITATION MECHANISMS AND DETECTABILITY 

We shall discuss three astrophysical processes that could be associated to the excitation 
of discontinuity g-modes and to the emission of gravitational waves: the coalescence of a 
NS-NS binary system, the onset of a phase transition in the neutron star core, and a pulsar 
glitch. For the last two cases, we shall estimate how much energy should be deposited into 
the modes in order the emitted signal to be detectable by interferometric detectors very 
sensitive in the kHz frequency region. 

The tidal excitation of (non rotating) NS oscillations in binary systems has been studied 
extensively in the past. In the case of core g-modes, it has been shown that the overlap 
integral between the displacement field and the tidal force field is small ( [Reisenegger fc| 
|Goldreich 1994| ; |Ho fc Lai 1999j ; [Lai 1994| ), and since the energy absorbed by a mode is 



proportional to the square of this integral, they are basically ineffective. To ascertain whether 
g-modes due to density discontinuities could be more efficiently excited in a binary system, 
we adopt the formalism introduced in previous works QGualtieri et al. 200 1| ; [F ons et al 



2002a| ), where we computed the gravitational perturbations of a neutron star induced by an 
orbiting companion by using a perturbative approach. We find that the typical resonance 
width of a discontinuity g-mode (Az/ fh 10~ 2 Hz) is much smaller than that of the f-mode 
(Av ~ 100 — 200 Hz). Similarly to what happens for the g-modes arising from thermal 
and composition gradients, the shift that the resonant excitation of a discontinuity g-mode 
produces in the number of cycles of the emitted signal is of the order of 10~ 2 cycles, too 
small to affect the detection significantly. For this reason, tidal excitation of g-modes during 
the coalescence of NS binary systems does not seem to be an efficient mechanism for the 
emission of gravitational radiation, and we shall not consider it in the following. 

A second excitation mechanism is the onset of a phase transition. According to current 
theories of stellar structure and evolution, the matter composing a neutron star may undergo 
a phase transition to quark matter if, at some point of the star life, a sufficiently high density 
is reached in its core. The critical density required to produce such a transition can be reached 
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soon after the core-collapse which gives birth to the neutron star, or at a later stage of the 
evolution, due to accretion ( Cheng fc Dai 1998| ) and/or spin-down ( [Ma fc Xie 1996|) . The 
onset of a phase transition is accompanied by a sudden reduction of the stellar radius and 
the difference in binding energy between the initial and final configuration is expected to 
be radiated away. The binding energy released during this late mini-collapse is of the same 
order as that emitted during neutron star formation following a gravitational core-collapse, 
i.e. « 2 x 10 53 erg = O.lM c 2 QAlcock, Farhi fc QlintoT986| ; |Bombaci fc Datta 2000|) . How 



much of this energy reservoir can be redistributed in the different modes of oscillation of the 
star and radiated in gravitational wave is, however, an open issue. Therefore, the only thing 
we can reasonably estimate is which fraction of this energy should go into the quasi normal 
modes in order the signal emitted at the corresponding frequencies to be detectable. 

The quasi normal modes of oscillation of a star may also be excited as a consequence of a 
glitch. Glitches are sudden changes in the rotation frequency of the neutron star crust which 
superimpose to the usual gradual spin-down under magnetic torque. They are observed in 
many pulsars and are thought to be related to quakes occurring in the solid structures such 
as the crust, the superfluid vortices and, perhaps, the lattice of quark matter in the stellar 
core ( |Anderson fc Itoh 1975| ; |Kuderman 1976| ; |Pines fc Alpar 1985| ). The observed glitches 
are very small, with a typical fractional spin variation of about AQ/Q pa 1CT 6 — 10~ 8 , which 
would be associated to an energy release of the order of AE ~ IQAQ, where I is the 
moment of inertia of the star. For the glitches observed in the Crab and Vela pulsars this 
simple estimate gives AE ~ 2 • 10~ 13 M c 2 and AE ~ 3 • 1CT 12 M Q c 2 , respectively. As 
before, we shall assume that a fraction of this energy is redistributed into the quasi normal 
modes of oscillation of the perturbed star, and we shall evaluate how much is needed to 
excite a mode to a detectable level. 

4.1 Detectability 

The frequency of the discontinuity g-modes we have studied ranges between u g ~ [0.5 — 
1.4] kHz, while Uf ~ [1.3 — 3.2] kHz and for the first p-mode v v ~ [2.5 — 7.3] kHz. In order 
to detect a signal emitted as a consequence of the excitation of these modes a very sensitive 
high frequency detector would be needed. In the following we shall consider the gravitational 
laser interferometric detector EURO, that has recently been proposed in a preliminary as- 
sessment study and for which the sensitivity has been estimated by Sathyaprakash and 
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Schutz ( [littp : / / www. astro . cf. ac.uk/ geo / euro|) . EURO should be extremely sensitive in the 



kHz region, which is crucial to infer informations about neutron stars interior from gravi- 
tational wave data, and it has been envisaged in two experimental configurations: one for 
which the noise curve in the frequency region 10 — 10 4 Hz can be written as 



SJu) = 10 



-50 



3.6 10 9 1.3 10 5 o / u 2 ^ 

-^- + -^- + 1.3 1(T^ + 



Hz" 1 , (16) 

V* V \ If ' 

where = 10 3 Hz; the second configuration is the more ambitious Xylophone detector, for 
which the shot noise, i.e. the third term within brackets in (|16"D, should be eliminated by 
running several narrow-banded interferometers. We will refer to these two possible configu- 
rations as EURO and EURO-Xylo, respectively. We shall assume that the gravitation signal 
from a star pulsating in a quasi- normal mode is a damped sinusoid of the form ( Echeverria 



h(t) = Ae {t ^- t)/r sin [2tiu (t - t m )] , (17) 

where t arr is the arrival time and r the damping time of the oscillation. The amplitude A 
can be expressed in terms of the total energy AE Q = AE/M Q c 2 deposited into the mode, 
and of the distance D of the source by ([Kokkotas, Apostolatos fc Andersson 2001|) 



V 



Given a detector with noise spectrum S n (v), and assuming that the matched filtering tech- 
nique is used to extract the signal, the signal to noise ratio (SNR), can be written as 

r a 2 tV /2 

5M= f— , (19) 

where T = 4Q 2 /(1 + AQ 2 ) is a form factor, and Q = ttut is the quality factor of the 
oscillation. 

It should be noted that in the case of discontinuity g-modes, since the gravitational 
damping time ranges between r w (10 7 — 10 12 ) s, the oscillations will be damped by other 
dominant dissipative mechanisms rather than by the emission of gravitational waves. In 
newly born neutron stars (or hybrid stars with a quark core) the dominant damping mech- 
anism is neutrino viscosity, with a typical timescale of 10 — 10 2 s ( |van den Horn fc van 



[Weert 1981| ; Goodwin fc Pethick 1982] ; [Thompson fc Duncan 199~3| ). In the glitches scenario 



(or other similar excitations occurring in late stages of the evolution of neutron stars) the 
neutron star is cold and neutrinos are not coupled to the matter. The dissipative mecha- 
nisms are then associated to the kinematical properties of matter. In this case, bulk, shear 
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viscosity and boundary layer dissipation are more efficient than gravitational damping. In 
any case, the typical dissipative time-scales are thought to be much larger than the neu- 
trino damping time-scale in proto-neutron(quark) stars ( [Bildsten fc Ushomirsky 2000| ), but 
a word of caution is needed. It has been found that quark-matter viscosity is much larger 
than that of normal neutron star matter ( |Madsen 1992| ), and therefore it would quickly 
damp the oscillations in the dense quark core (r ~ 1CT 2 s). However, these results are still 
under debate, and it remains to be carefully analyzed how the outer shell would react to the 
strong damping in the inner region, and whether or not boundary effects would change the 
picture. Note however that, for the range of frequencies considered in this paper, the quality 
factor Q ^> 1 as long as r is greater that, say, 10~ 2 s; thus, the form factor in Eq. ([L9|) can 
be taken equal to one, and the signal to noise ratio happens to be quite insensitive to the 
damping time. 

In the following, we shall estimate the amount of energy that should go into a quasi 
normal mode to produce a gravitational signal detectable with a signal to noise ratio SNR = 
3, either by EURO or by EURO-Xylo detector. For g-mode pulsations, we consider the two 
limiting frequencies V\ = 0.5 kHz and v% = 1.4 kHz, being the g-mode frequency of all our 
models included in this range. The excitation of a core g-mode at 0.1 kHz ( [Reisenegger & 
Goldreich 1992| ; |Lai 1994]) is also computed for comparison. For the f- and p-modes, our 



estimates are done for two frequencies, V[ = 2.3 kHz and u v = 5.5 kHz, which are obtained 
as an average over the corresponding range of variation for all our models. The gravitational 
damping time of the fundamental and p-mode are smaller than the other dissipative time 
scales; thus, in this case the emission of gravitational waves is the dominant dissipative 
mechanism and both T{ and r p have been used explicitly in our estimates. 

For the two different excitation scenarios, mini-collapse and glitches, we shall consider 
sources located at distances D = 15 Mpc and D — 10 kpc, respectively. Note that the 
Virgo cluster, with its some 2 10 3 members, is at D m 15 Mpc and that the Vela and Crab 
pulsars, that are known to produce glitches, are at distances of D = 0.5 kpc and D = 2 kpc, 
respectively. 

Our estimates are given in Table (|5J), and show that the high frequency detectors would be 
a very interesting instrument to probe, through gravitational signals, the nature of neutron 
stars interiors. Note that, since for the EURO configuration the sensitivity at high frequency 
is strongly limited by the shot noise, a p-mode at u p = 5.5 kHz could be detected only with 
the advanced EURO-Xylo configuration. If one focuses on the fundamental mode and the 
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discontinuity g-mode, the minimal energies required to detect a signal with SNR = 3 with 
the EURO detector range within 1.6 ■ 1(T 13 < AE/M Q c 2 < 3.0 • 10" 5 , depending on the 
source distance and on the frequency of the signal. With the EURO-Xylo configuration, 
this range becomes 4.1 • 10~ 14 < AE/M Q c 2 < 3.4 • 1CT 7 . As mentioned before, the amount 
of binding energy which is expected to be released in a mini-collapse is approximately the 
same as that emitted during a core-collapse Supernova. In the Supernova case, it has been 
suggested that the fraction of this energy radiated in gravitational waves is in the range 
[1(T 8 - 1(T 6 ] M c 2 QMonchmeyer et al. 199H Pimmelmeier, Font fc Miiller 200l| ; [Zwerger 



|fc Miillcr 19971 ); if we assume that a comparable amount of energy goes in gravitational 
waves also in the case of a mini-collapse, the chance for a quasi normal mode signal to be 
detectable from a source closer than 15 Mpc is promising. The situation is clearly even more 
interesting for galactic sources, since the amount of energy required to excite a mode at a 
detectable level is typically a small fraction of the total energy reservoir. 

5 CONCLUSIONS 

According to the relativistic theory of stellar perturbations, a star perturbed by any external 
or internal process emits gravitational waves at the characteristic frequencies of its quasi 
normal modes. It has been shown that both the mass and the radius of a pulsating cold 
neutron star could be inferred, if the quasi normal modes were identified by a spectral 
analysis in the detected wave flAndcrsson fc Kokkotas 1998|) . More recently, the possibility 
of tracing the presence of a superfluid core through the detection of superfluid modes has 
also been investigated QAndersson fc Comer 2001" ). As described in the introduction, g- 



modes are also good markers of the internal composition of a neutron star: low frequency 
g-modes indicate a non homogeneous composition in the outer layers or in the core of the 
star, or a thermal profile, whereas the high frequency g-modes which we study in this paper, 
are associated to phase transitions occurring at supranuclear densities (Sotani et al, 2001). 
Although the analysis has been restricted to simplified polytropic EOSs, our results can be 
generalized to realistic equations of state, since it is known that the oscillation frequencies 
depend more on global properties, like mass and radius, than on the specific form of the 
microscopical interactions. We have analyzed models of stars with a mass of 1.4 M Q , since 
astrophysical observations show that most of the well measured masses of neutron stars in 
binary systems cluster in a narrow range around this value. 
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We find that high frequency discontinuity g-modes exhibit several interesting features: i) 
the linear relation known to exist between the frequency of the fundamental mode and the 
square root of the average density, still holds for neutron stars with density discontinuities. 
Thus, a measure of the f-mode frequency would be useful to constrain the neutron star 
radius, provided the mass is measured independently, even when a phase transition takes 
place in its core. However, the presence of the discontinuity introduces a larger error on 
the NS radius determination, ii) Should density discontinuities actually exist, the associated 
g-mode will have a frequency of about 1 kHz, somewhat lower than that of the f-mode, 
and therefore clearly distinguishable. It may be reminded that thermal or chemical g-modes 
appear at significantly lower frequencies (< 200 Hz), iii) Since the g-mode frequency strongly 
depends on the amplitude of the discontinuity, it could be used to infer the value of Ap/p d . 
In addition, since the maximum of v g , considered as a function of the discontinuity radius 
i?d, is related to Ap/p<j and to the average density, p i; of the core internal to Rd, (see eqs. 
(IH) and (pT5[)) it would be possible to set a lower bound on pi, and to rule out all EOSs for 
which the combination C 2 = p\ is smaller than the observed one. 

In the spirit of previous works ( [Andersson fc Comer 20TJT| ) , which discussed the ob- 
servability of oscillation modes produced by the existence of superfluid components in the 
interior of a neutron star, we have analyzed some astrophysical processes in which high 
frequency modes could be excited, and we have evaluated the amount of energy that the 
modes should convey to produce a gravitational signal detectable with a SNR=3 by the high 
frequency gravitational detectors under study. We find that the required energy does not 
conflict with the estimates of the energy expected to be released in gravitational waves at 
the onset of a phase transition, or during a glitch. In particular, the detector EURO, espe- 
cially in the configuration known as EURO-Xylophone, would have a significant chance to 
reveal the gravitational signal associated to the excitation of the fundamental mode, of the 
first pressure mode and of the discontinuity g-mode following a galactic glitch or a stellar 
core-collapse or a mini-collapse induced by the onset of a phase transition to quark matter 
within 10-15 Mpc. By cross-correlating the information carried by these mode frequencies, 
and in particular by the high frequency discontinuity g-modes, we would acquire extremely 
valuable indications on the properties of nuclear matter at supranuclear density. 
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R (km) 


p c (g/cm 3 ) 


\fp (km r ) 


Pd (g/cm 3 ) 




u { (kHz) 


u p (kHz) 


v s (kHz) 


13.44 


0.92 10 15 


0.0292 




0.0 


1.666 


4.045 




12.04 


1.39 10 15 


0.0344 


3 10 14 


0.1 


1.998 


4.637 


0.504 


12.22 


1.36 10 15 


0.0337 


4 10 14 


0.1 


1.962 


4.548 


0.567 


12.42 


1.32 10 15 


0.0329 


5 10 14 


0.1 


1.915 


4.459 


0.613 


12.65 


1.27 10 15 


0.0319 


6 10 14 


0.1 


1.857 


4.365 


0.644 


12.92 


1.19 10 15 


0.0310 


7 10 14 


0.1 


1.792 


4.260 


0.659 


13.21 


1.11 10 15 


0.0300 


8 10 14 


0.1 


1.723 


4.146 


0.658 


13.43 


1.02 10 15 


0.0292 


9 10 14 


0.1 


1.670 


4.052 


0.641 


10.71 


2.17 10 15 


0.0410 


4 10 14 


0.2 


2.408 


5.325 


0.840 


10.99 


2.07 10 15 


0.0395 


5 10 14 


0.2 


2.330 


5.155 


0.912 


11.35 


1.95 10 15 


0.0376 


6 10 14 


0.2 


2.226 


4.970 


0.961 


11.83 


1.77 10 15 


0.0354 


7 10 14 


0.2 


2.088 


4.750 


0.987 


12.50 


1.51 10 15 


0.0325 


8 10 14 


0.2 


1.901 


4.452 


0.979 


13.38 


1.15 10 15 


0.0294 


9 10 14 


0.2 


1.680 


4.072 


0.906 


8.68 


4.44 10 15 


0.0562 


5 10 14 


0.3 


3.216 


6.683 


1.211 


9.13 


3.95 10 15 


0.0521 


6 10 14 


0.3 


3.039 


6.325 


1.281 


9.66 


3.46 10 15 


0.0479 


7 10 14 


0.3 


2.831 


5.953 


1.326 


10.41 


2.88 10 15 


0.0428 


8 10 14 


0.3 


2.553 


5.506 


1.339 


12.15 


1.87 10 15 


0.0340 


9 10 14 


0.3 


2.002 


4.632 


1.251 



Table 1. The frequencies of the fundamental mode, of the first p— mode and of the discontinuity g-mode are shown for a set 
of stellar models with M = 1.4 Mq and T = 2. For each star, the equilibrium parameters are also tabulated. The polytropic 
coefficient K is K(l + Ap/p^) 2 = 180 km 2 , so that all stars have the same low-density equation of state. The first entry 
corresponds to a model with no density discontinuity. 
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R (km) 


p c (g/cm 3 ) 


yfp (km r ) 


Pd (g/cm 3 ) 


Ap/Pd 


Vf (kHz) 


Up (kHz) 


v s (kHz) 


17.96 


0.72 10 15 


0.0189 




0.0 


1.286 


2.543 




13.38 


2.14 10 15 


0.0294 


4.0 10 14 


0.1 


2.035 


3.747 


0.641 


14.05 


1.86 10 15 


0.0273 


4.5 10 14 


0.1 


1.895 


3.528 


0.648 


14.75 


1.61 10 15 


0.0254 


5.0 10 14 


0.1 


1.764 


3.321 


0.650 


15.51 


1.38 10 15 


0.0235 


5.5 10 14 


0.1 


1.633 


3.311 


0.646 


16.39 


1.15 10 15 


0.0217 


6.0 10 14 


0.1 


1.495 


2.892 


0.635 


17.32 


0.94 10 15 


0.0199 


6.5 10 14 


0.1 


1.367 


2.680 


0.614 


17.90 


0.81 10 15 


0.0190 


7.0 10 14 


0.1 


1.293 


2.555 


0.593 



Table 2. The equilibrium parameter and the quasi normal mode frequencies are shown as in Table 1 for stellar models with 
M = 1.4 Mq and T = 1.67. The polytropic coefficient K is chosen in such a way that the star with no discontinuity (Ap/pd = 
) and those with Ap/pd = 0.1 have the same low-density equation of state. 



R (km) 


Pc (g/cm 3 ) 


yfp (km 1 ) 


Pd (g/cm 3 ) 


Ap/Pd 


!/ f (kHz) 


f p (kHz) 


v g (kHz) 


15.17 


0.82 10 15 


0.0243 




0.0 


1.490 


3.308 




10.22 


3.40 10 15 


0.0440 


4.0 10 14 


0.2 


2.736 


5.361 


0.932 


10.51 


3.16 10 15 


0.0422 


4.5 10 14 


0.2 


2.639 


5.185 


0.961 


10.82 


2.94 10 15 


0.0404 


5.0 10 14 


0.2 


2.542 


5.016 


0.985 


11.17 


2.72 10 15 


0.0385 


5.5 10 14 


0.2 


2.434 


4.839 


1.003 


11.58 


2.48 10 15 


0.0365 


6.0 10 14 


0.2 


2.314 


4.648 


1.014 


12.07 


2.22 10 15 


0.0343 


6.5 10 14 


0.2 


2.177 


4.434 


1.016 


12.71 


1.92 10 15 


0.0317 


7.0 10 14 


0.2 


2.011 


4.174 


1.007 


13.66 


1.53 10 15 


0.0285 


7.5 10 14 


0.2 


1.789 


3.817 


0.972 


15.08 


1.03 10 15 


0.0246 


8.0 10 14 


0.2 


1.507 


3.338 


0.874 



Table 3. The equilibrium parameter and the quasi normal mode frequencies are shown for stellar models with M = 1.4 Mq 
and T = 1.83. The first entry is a star with no discontinuity, the other stars have a discontinuity of amplitude Ap/pd = 0.2, 
and the polytropic coefficient K is chosen as in Table 2. 
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R (km) 


p c (g/cm 3 ) 


yfp (km x ) 


Pa (g/cm 3 ) 


Ap/Pd 


u { (kHz) 


Up (kHz) 


v s (kHz) 


11.32 


1.22 10 15 


0.0378 




0.0 


2.005 


5.361 




8.53 


3.44 io 15 


0.0577 


6.0 10 14 


0.3 


3.095 


7.329 


1.171 


8.71 


3.34 10 15 


0560 


7.0 10 14 


0.3 


3.031 


7.145 


1.253 


8.92 


3.21 10 15 


0.0540 


8.0 10 14 


0.3 


2.945 


6.949 


1.317 


9.20 


3.03 10 15 


0.0516 


9.0 10 14 


0.3 


2.831 


6.734 


1.362 


9.58 


2.77 10 15 


0.0485 


1.0 10 15 


0.3 


2.673 


6.473 


1.395 


10.14 


2.41 10 15 


0.0445 


1.1 10 15 


0.3 


2.445 


6.116 


1.368 


10.60 


2.12 10 15 


0.0417 


1.15 10 15 


0.3 


2.268 


5.829 


1.326 


11.26 


1.66 10 15 


0.0381 


1.2 10 15 


0.3 


2.027 


5.404 


1.208 



Table 4. The equilibrium parameter and the quasi normal mode frequencies are shown as in previous Tables, for stellar models 
with M = 1.4 Mq, T = 2.25 and Ap/p d = 0.3. 



mode i^modc D = 15 Mpc D = 10 kpc 

AE (Euro) AE (Euro-Xylo) AE (Euro) AE (Euro-Xylo) 



gcore 


0.1 kHz 


3.5 


io- 


-7 


M Q c 2 


3.4 


10" 


-7 


Mqc 2 


1.6 


10" 


13 


Mqc 2 


1.5 


10" 


13 


Mqc 2 


Sdisc 


0.5 kHz 


3.9 


io- 


-7 


M q c 2 


1.0 


10" 


-7 


Mqc 2 


1.7 


10" 


13 


Mqc 2 


4.5 


10" 


14 


MqC 2 




1.4 kHz 


5.4 


10" 


-6 


M @ c 2 


9.2 


10" 


-8 


Mqc 2 


2.4 


10" 


-12 


Mqc 2 


4.1 


10" 


14 


Mqc 2 


f 


2.3 kHz 


3.0 


10" 


-5 


M Q c 2 


9.2 


10" 


8 


Mqc 2 


1.4 


10" 


11 


Mqc 2 


4.1 


10" 


14 


Mqc 2 


pi 


5.5 kHz 


8.7 


io- 


4 


Mqc 2 


9.1 


io- 


8 


Mqc 2 


3.8 


10" 


10 


Mqc 2 


4.0 


10" 


14 


Mqc 2 



Table 5. The amount of energy that should be radiated in gravitational waves at the frequency of a given mode, in order 
for the signal to be detectable with a SNR = 3 by the EURO detector configurations. The source distance is taken to be 
D = 15 Mpc and D = 10 kpc (see discussion in the text). 
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2 




q 2 < 1 1 1 1 1 1 1 1 1 - 

8 9 10 11 12 13 14 15 16 17 

R ( km ) 

Figure 1. Wc compare the mass-radius curve of stellar models without discontinuity (continuous line) and with a high density 
discontinuity in the core (dashed line). Both curves refer to the same low-density equation of state with T = 2, the only 
difference being the presence of a discontinuity at p<j = 7 • 10 14 g/cm 3 with Ap/p^ = 0.3, for the model shown with a dashed 
line. As discussed in the text, when a discontinuity is present the star is more compact and the maximum mass is lower. 
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Figure 2. The f-mode and the first p-mode frequencies are plotted as a function of the density at the discontinuity, p^, for 
the models of stars with mass M = 1.4 Mq and T = 2 given in Table |lj The different curves refer to assigned values of 
Ap/Pd = ( 0.1, 0.2,0.3 ), increasing from bottom to top. The variation of the mode frequency with the discontinuity parameters 
depends on the associated variation of the mean density and the compactness of the star (see text). All mode frequencies tend 
to the same limiting value (black dot) as approaches the central density, which is the frequency of the mode when there is 
no discontinuity. 
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Figure 3. We show the f-mode frequency as a function of the mean density of the NS for a sample of our models, labelled 
with (r, Ap/p d ). The fit by Andersson and Kokkotas (1998), obtained for stars without density discontinuities, is shown as a 
continuous line (AK98-fit), while the black dots are the f-mode frequency of our models with Ap/pd = 0. The Figure shows 
that, for each assigned value of Ap/pd, the f-mode frequency is still a linear function of \J M/R s , as in the case of stars with 
no discontinuity, but with a different slope. 
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Figure 4. The g-mode frequency of the stellar models given in Tables is shown as a function of Pd- It is remarkable that 
u g depends much more on the amplitude of the density discontinuity Ap/pd, than on the other parameters of the EOS. Should 
a g-mode ever be detected, this feature would allow to infer the value of the density jump. 
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Figure 5. The discontinuity g-mode frequency is plotted as a function of the normalized radius of the internal core for two 
selected stellar models with T = 1.83 and T = 2, respectively. The amplitude of the density jump is Ap/p d = 0.2. The maximum 
frequency is reached when the radii of the internal core and of the outer region are comparable. For all the studied models 
I'max IS reached for 0.4 < R d /R < 0.6. 
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